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Abstract 

Using a modified perturbative approach that includes the Sudakov resummation 
and transverse degrees of freedom we analyze the penguin-induced B~ — > 7r~</> 
decay by applying the next-to-leading order effective weak Hamiltonian. The mod- 
ified perturbative method enables us to include nonfactorizable contributions and 
to control virtual momenta appearing in the process. Besides, we apply the three- 
scale factorization theorem for nonleptonic processes that offers the possibility of 
having the scale-independent product of short- and long-distance parts in the am- 
plitude of the weak Hamiltonian. The calculation supports the results obtained in 
the BSW factorization approach, illustrating the electroweak penguin dominance 
and the branching ratio of order O(10~*). However, the estimated prediction of 
16% for the CP asymmetry is much larger than that obtained in the factorization 
approach. 
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I. INTRODUCTION 



Among a variety of heavy-meson decaying channels, exclusive two-body nonleptonic decays are 
theoretically the most challenging ones, owing to the phenomenon of hadronization and the effects 
of final-state interactions. On the other hand, they present the most promising way to detect CP 
violation in the heavy-meson sector and to explore the CKM mixing matrix elements. 

The mechanism of CP violation can be investigated directly in the charged sector of B mesons 
by measuring CP asymmetry. CP asymmetry is defined as a relative difference between of the 
decay rates of the B meson and its CP-conjugated state, i.e., 

_ r(B- ^ /) - T(B- ^ 7) 

r(B-^/) + r(B+^/) ^ ' 

Nonvanishing CP asymmetries appear through the interference between amplitudes with different 
weak CP-violating phases and different CP-conserving strong phases coming from the final-state 
strong interactions different from zero. 

Nowadays, experimental facilities offer a possibility of searching for CP asymmetries in penguin- 
induced nonleptonic decays, the very promising decays to detect direct CP violation. Such decays 
have small branching ratios (BR), but satisfy both requirements for CP- violating asymmetry, owing 
to the fact that penguins are loop diagrams with different quark generations contributing with dif- 
ferent weak CP-phases from the CKM matrix and that final-state strong interaction phases emerge 
from the absorptive part of penguin amplitudes. This mechanism of generating CP asymmetries 
in decays that involve penguins was first considered by Bander, Silverman, and Sony 

In this paper we discuss the pure penguin- induced B~ n~(f) decay governed by the heavy- 
quark b dss decay. Performing a consistent 1/Nc expansion {Nc, quark color number) of QCD- 
penguin amplitudes, one can show that in this process the QCD-penguin contribution should be 
suppressed and the dominant contribution comes from the electroweak (EW)-penguin operators 
0. The B^ — > 71^ (f) process has already been considered by many authors [^- within the Bauer- 
Stech-Wirbel (BSW) factorization approach [0. This method for reducing the hadronic matrix 
element of four-quark operators to the product of two current-matrix elements cannot account 
for QCD interactions between the currents, except by parametrizing them by a phenomenological 
parameter in the generalized factorization approach In general, momenta of virtual gluons or 
photons in a process, appearing explicitly in the penguin matrix elements after factorization have 
to be considered free parameters. CP asymmetry depends strongly on these parameters and the 
predictive power of calculations performed within the factorization prescription is greatly reduced. 

Our aim is to investigate the B~ — 7r~0 decay in the modified perturbative approach. Per- 
turbative calculations of exclusive B decays were carried out by different authors all of them 
following the framework for analyzing exclusive decays in the perturbative QCD (pQCD) ap- 
proach developed by Brodsky and Lepage, and other authors [|1^]. In the perturbative approach, 
exclusive amplitudes involving large momentum transfer factorize into a convolution of a process- 
independent and perturbatively incalculable distribution amplitudes (hadronic wave functions), 
one for each hadron involved into the decay, with a process-dependent and perturbatively calcula- 
ble hard scattering amplitude of valence partons. 

The apphcability of such a pQCD framework to exclusive decays was widely discussed [|TT|, 



12] owing to the concern about the possible uncontrollable nonperturbative (end-piont region) 



contributions and the problem was solved in a modified perturbative approach proposed by Li and 



Sterman ITS 
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Besides offering a reliable perturbative calculation, the modified perturbative approach offers a 
possibility of going beyond the factorization approximation in the calculation of four-quark matrix 
elements. It also enables us to assign the pro cess- dependent virtual momenta in the loop 
matrix elements and to fold them with their distribution in a particular decay. In this way, the 
uncertainties in CP asymmetry, which are due to some ad hoc quark model values of as applied 
in the factorization approach, do not appear. 

The purpose of this paper is to present a complete calculation of factorizable and nonfactor- 
izable contributions in the penguin-induced B~ iT~(j) decay up to order O(asaem), testing the 
results on various CKM mixing matrix parameters. Especially, we wish to examine nonfactor- 
izable contributions from the QCD-penguin operators and assign their role in the EW-penguin 
dominated processes, such as B~ 7r~0. 

The plan of the paper is as follows. In Sec. II we introduce the method of calculation based 
on the next-to-leading order (NLO) effective weak Hamiltonian, and the modified perturbative 
approach. A detailed analysis of the 7r^0 process is presented in Sec. III. A discussion of 

mesonic wave functions and the Sudakov form factors is given in Sec. IV and the selection of proper 
mesonic wave functions is made. In Sec. V we present our numerical results for the branching 
ratio and CP asymmetry, comparing them with those obtained from the factorization approach 
and examine the dependence of the results on the choice of CKM parameters. Concluding remarks 
are given in Sec. VI. 



II. PERTURBATIVE MODEL FOR CALCULATING THE 5 ^ vr DECAY 

The nonleptonic B~ — > 7i~(f) decay is governed by the weak decay of the heavy b-quark, b dss. 
The light antiquark of the B meson is the spectator in the decay, being only slightly accelerated 
by the exchange of a hard gluon to form a pion in the final state. 

In this section we present the basic ingredients for a calculation of such a penguin- induced 
decay. The first ingredient is the NLO effective weak Hamiltonian, which allows a consistent study 
of nonleptonic decays in which penguin operators are involved. The second ingredient, on which 
the paper is based, is a modified perturbative method for calculating exclusive decays by which 
the matrix elements of four-quark weak-Hamiltonian operators are perturbatively calculable. 



A. Low-energy effective weak Hamiltonian beyond the leading logarithmic approximation 



Following Ref. [14| we consider the NLO effective weak Hamiltonian for b ^ d transitions: 



H^si^B = -1) = ^ ^ vj c,ifi)0? + C2ifi)Oi'^ + E cMO, . (2.1) 

V 2 q=u,c \ k=3 ) 

The scale-dependent Wilson coefficients Cj(yu) are the short- distance part of the Hamiltonian 
and include NLO QCD corrections and leading-order ctem corrections. With Vq we denote products 
of CKM mixing matrix elements relevant to 6 — >■ c? transitions, Vq = V*^Vqb. Local four-quark 
operators, renormalized at the scale n, are 

0[''^ = {daqp)v-A{q(3ba)v-A , 02'' = {dq)v~A{qb)v-A , 

O3 = (db)v-A J2(^l')v-A , = {dahl3)v-A E(?/3?'a) V-A , 



3 



b 



d 



FIG. 1. QCD and QED one-loop penguinlike contributions of the tree-level operators Oi and O2 in 
the h — > dq'q' decay. 

C5 = (db)v-A E(W)v+^ , Oe = (dabf3)v~A E(?/3?'a) v+a , 

q' q' 

3 — — 3 - — 

q' q' 

3 — — 3 - — 

C9 = ^idb)v-Aj2^g'(^'^')v-A , Oio = -{dabp)v-Aj2^i'(9'l39'c,)v-A , (2.2) 

q' q' 



where V ± A = 1/2 7^(1 ± 75), a and P are color indices, q' G {m, (i, s, c, 6}, and e,/ are the 
corresponding quark charges. 0[^^ and 02''^ are tree-level operators, 0^,...,0q are QCD-penguin 
operators, and O7, Oio are EW penguins. 

From the quark content of the operators O^f* and 0^2^ it is obvious that they do not contribute 
at the tree level in 6 — dss transitions. Such transitions are pure penguin-induced, receiving 
contributions from the operators O3, Cio, in which q' is restricted to be a strange quark, q' = s. 

In the NLO weak Hamiltonian ( |2.1| ) the renormalization-scheme dependence of the Wilson 
coefficients is explicitly canceled by the inclusion of the one-loop QCD and QED matrix elements 
of the tree-level operators 0[^l, Fig. |l|. The one- loop matrix elements of the NLO weak Hamiltonian 
for a 6 — > dss transition can be written in terms of products of the tree-level matrix elements of 
penguin operators and the renormalization scheme-independent coefficients c(/i) : 

{dss\n.ff{AB = -i)\b) = ^ E VqiY^cMW- 

y Z q=u,c \k=3 

+^(3ci(/^) + C2(/i)) {((y - AGK,g^^2)) (O, + 0,)y^\] , (2.3) 



where (0^)*"''=^ = {dss\Ok\b) 



tree 



The function AG{m?, g^, /x^) arises from the one-loop penguinlike diagrams with q = u,c quarks 



in the loop: 

"1 



AG(m^,g^/x^) = -4 duu{l - u)\n ""M . (2.4) 
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Here, as well as the quark mass ruq and the renormalization scale /i, there appears a new parameter 
g^, which is the momentum squared of a virtual gluon or photon emerging from the loop. In the 
b dss transition can be indentified with the sum of the strange quark momenta squared. 

Concentrating now on the specific process B^{bu) —>■ TT~{du)(f){ss), a few comments are in order. 
The strange quark and antiquark building 0-meson are obviously in the color-singlet state. The 
ss pair coming from the virtual gluon decay in the QCD-penguin diagrams builds a color-octet 
state. Therefore, first, one expects a small contribution from the QCD-penguin operators and, 
second, the one-loop QCD penguinlike contribution (shown in Fig. |l| with an exchanged gluon) is 
not present in the 7r^0 decay. The final expression for the matrix element is then 

{n-meffiAB = -i)\B-) = ^ E vAy^cmw- 

y I q=u,c \k=3 

+ ^ - AGK, (g^),^)) (3ci(^) + C2(/i)) {(Or)*^^^ + (Os)*^^^}) • (2.5) 

Here {Ok)'''''' ~ {7r-(j)\Ok\B-) . 

We have retained the ttem-pi'opotional term as a part of the NLO weak Hamiltonian, although 
we show later that, in the perturbative approach to the order we are working with, such a contribu- 
tion emerges naturally. This term is producing hard final-state interaction phase shifts, necessary 
for generating CP asymmetry, which are due to the on-shell quarks rescattering in the loop for 
particular values of q^. The CP asymmetry depends strongly on the value of g^ and, in the 
factorization approach, the major source of uncertainties in predictions comes from the lack of 
information about the g^ value after the factorization of hadronic matrix elements is performed. 
On the contrary, in the perturbative approach, the g^-dependence of the loop amplitude is calcu- 
lable directly as a part of the hadronic matrix element and it is determined by the momentum 
distributions in a particular process. However, in the strict factorization, and in the perturbative 
calculation for the process considered, the average g^-value can be simply determined to be the 
mass of the cf) meson squared, (g^) = M^. Further discussion about this point is left for Sec. V. 

Let us now continue with the estimation of the matrix elements of four-quark operators. 

B. Modified perturbative approach to the calculation of the matrix elements of four-quark 

operators 

Perturbative calculations of matrix elements in exclusive hadron decays can be carried out in 
the Brodsky-Lapage (BL) formalism [|l^ . Hadrons are considered in the leading approximation as 
a bound state of valence quarks and/or antiquarks, depending on the hadron. The amplitude of 
the process factorizes into the convolution of distribution amplitudes of hadrons involved in the 
decay (hadron wave functions) and the hard scattering amplitude of valence partons. Hadronic 
wave functions represent the nonperturbative part, which has to be determined for heavy hadrons 
in relativistic constituent or nonrelativistic models, or for light hadrons by the QCD sum-rule 
method or by lattice calculations. Further discussion about the wave functions is given in Sec. IV. 

The hard scattering amplitude can be calculated perturbatively, taking into account all possible 
exchanges of a hard gluon between valence partons in a given a<j-order of the calculation. Valence 
quarks are carrying some fraction of momentum of their parent hadron and the final expression is 
integrated over the fractions Xi {i = 1, 2, 3), see Fig. |^. The concern about the applicability of this 
method to exclusive processes was raised when it was noted that even at large momentum transfer 
the contribution to these processes could come predominantly from the momentum regions in which 
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Xipi+ki X3P3+k3 




FIG. 2. The basic graph of the B tt (j) decay with the momentum definitions specified. The black 
circle stands for the NLO effective weak Hamiltonian 



Us is large |]TT|, ||T2|. The problem is connected with the end-point region of momentum fractions. 
Namely, when one of the hadron constituents carries all the momentum of a hadron, the situation is 
no more perturbative and significant uncontrollable soft contributions might appear. The solution 
of the problem was proposed by Li and Sterman |^ . Contrary to the BL formalism, they suggested 
to go beyond the collinear approximation and to retain the small transverse momentum of valence 
quarks. Owing to the transverse degrees of freedom, the parton virtualities become large enough in 
the whole region for a reliable perturbative calculation. Furthermore, they included the Sudakov 
form factor for each of the hadrons in the decay to suppress the contributions from dangerous 
soft regions. All these effects can be incorporated into the factorization formula by expressing the 
transverse momentum variables in the Fourier transformed 6-space. 

The very last formula, which is used throughout this paper to calculate the matrix elements of 
four-quark operators from (p.5|) relevant to the B~ n^cj) decay, is 



{n-<P\Ok\B-)= [dx] 



fb 

4:71 



^';(X3, h)mx2, h) m{x}, {b}, MB)^ij(xi, 6i)e 



-S[{x},{b},MB 



(2.6) 



where Xj {i = 1,2,3) are fractions of longitudinal momenta of B, 0, and vr mesons, respectively. 
Analogously, bi denote the Fourier- transformed transverse momenta of these mesons, [dx] = 
dxidx2dx3 and {x} denotes the set of variables {xi, X2, x^}. Similarly, for b variables. 
and are the wave functions of the outgoing vr and mesons, and the decaying B meson, 
respectively. Tk is the hard scattering amplitude describing the B ^ Tr(f) decay at the one-loop 
level caused by one of the operators Ok ( |2.2| ). The exponential factor in formula ( |2.6|) is the 
Sudakov factor. Its explicit form is given in Sec. IV. 



III. CALCULATION OF THE B' vr ^ DECAY 

The basic graph representing the B^ rt^cf) decay in the perturbative approach is shown in 
Fig. |. 

In order to perform the calculation, we have to specify the hadronic momenta. The simplest 
choice can be made using the light-cone coordinates and taking the B meson to be at rest: 
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Pb=Pi = ^(1,1,0), pi = Ml, 
^0 = P2 = ^(l,r^O), pI = MI, 

P. = P3 = ^(0,l-r2,0), p^^O. (3.1) 

Here, Mb and are the masses of the B and the meson, respectively, and r is defined as their 
ratio, r = M^/Mb- The pion is taken to be massless. In addition to carrying some portion of the 
momentum of their parent meson, the valence quarks also carry some small transverse momentum 
k (see Fig. ^. 

Working in the leading order, the perturbative part, namely, the amplitudes of the four-quark 
operators ( p.Gj) , can be calculated from the Feynman graphs with all possible attachments of a hard 



gluon, shown in Fig. ^ As discussed later explicitly. Fig. ^(a) shows factorizable amplitudes of 
the decay. Figure ^(b) shows nonfactorizable amplitudes, not presented in the BSW-factorization 
approach. Contributions of one-loop induced penguinlike diagrams coming from the tree-level 
operators are shown in Fig. ^(c). These contributions can be taken into account immediately 
assuming the NLO weak Hamiltonian, but they emerge naturally in the pertubative model. Finally, 
Fig. ^(d) shows some additional diagrams that have to be included to perform a proper C(asaem) 
calculation. 

In full hadron wave functions one can split up spin wave functions from the rest; particularly, 

^ij(a;i, 6i) = -^( A + Mb)i5^<I>b{xi, h) , 

K{x„h) = -^i.h^Ki^sM) , (3.2) 

where He is the unit color matrix, is the polarization vector of the meson. The scalar wave 
functions $ are specific to each of the mesons and are discussed in Sec. IV. 

In this process it is enough to calculate QCD-penguin contributions, because EW-penguin 
operators are simply related to QCD penguins as 

- 2 

03 = {db)v-A(ss)v-A = —Og , 

36^ 

- 2 

04 = -{ds)v-Aisb)v~A = i^Oio , 

3e., 



O5 = idb)v-Aiss)v+A = ^C*? , 



2 

3e., 

2 



06 = -2{ds)sMsb)s-p = i^Os , (3.3) 

36^ 

where is the strange quark charge, = — l/3e. 

Applying the standard procedure for calculating Feynman diagrams, the hard scattering am- 
plitudes Tfc can be easily worked out by performing color and spin traces. They appear after the 
spin and color parts of the wave functions are added to the amplitude, leaving the scalar functions 
$ at the place of the full mesonic wave functions \E' in (|2.6| ). 
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a) 





b) 





"loop 



-*loop 





d) 



FIG. 3. Leading-order contributions to the B~ Tr~(j) decay: a) factorizable, b) nonfactorizable, 
and c) and d) penguinlike. The square stands for the penguin operators Ok, k = 3, ...,10, and the circle 
represents the tree-level operator 02- 
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All operators receive the same contribution from the factorizable diagrams, Fig. ^(a): 
Tfact = Ta + Tb , 

AT T 

(1 + X3(l - r^))e-p3 



Ta = -47ra,(/.)C;^^32M|r 



1 1 



Nr 



32Mir 



xi{l — r )e ■ pi ~ xie ■ ps 



DgD, 

1 1 



(3.4) 



where Cf is the color factor equal to Cf = Cp = 4/3 for the O3 (Og) and {O7) operators and 
Cf = Cp/Nc for the others. We are going to take the number of colors Nc to be equal to 3. 

Dg, Dh, and Dg^ denote the denominators of the exchanged gluon, and of the virtual b and d 
quark propagators, respectively: 



DG = ql + iv = Ml[xl - XiX3(l - r^)] - {k^ - ^3)^ + iy^ 
Db = ql + ir]- -Mlxsil - r^) - kl + ir] , 
Dq^ = ql +iv = Ml[xl - xi(l - r^)] -kl + ir] . 



(3.5) 



The approximation made here was to take m?, ~ Mb in the propagator of the heavy h quark. 
Other propagators are massless. 

The operators (Oio) and (Os) also receive nonfactorizable contributions coming from 
the estimation of the diagrams in Fig. ^(b). The expressions for the (Oio) operator are 



Tnonfact{Oi) — Tc{Oi) + T£,{Oi) , 

Tcm = -47r«,(/.)^^^32M|r 
Tnm = -47ra.(/x)^^^32M|r 
and similarly for the Oq (Og) operator 

TnonfactiOo) = Tc{Oq) + Td{Oq) , 



[1 - r^)(l - xi - X2)e-pi 



DgD,, 



(2xi - X2 - X3 - r^(x2 - X3))e ■ ^3 



1 1 



(3.6) 



_47ra,(/i)^^^32M|r 



3 73 2^2""^'^'' 
1 -2xi - X2 + X'i + r^(l - X2- X3))e ■ ps 



1 



Dg D, 



92 



[1 - r^)(xi - X2)e-pi 



3 V32v^""'' 

The denominators of the virtual quark propagators in Fig. ^b) are 



1 



Dg D 



(3.7) 



D 
D 



92 



93 



qq2 +iv = Ml (1 - xi - X3)(-(xi - X3) + {1~X2- X3)r'^) - {ki + k2- ksf + irj 
Qq. +iri = Ml {xi - X2){{xi - X3) - {x2 - x^y"^) - {ki - h - hf + irj . 



(3. 



The calculation of the one-loop EW-penguinlike contributions from Fig. ^(c) follows the already 
familiar procedure. Both diagrams Aioop and -Bjoop are propotional to their skeleton graphs A and 
B , respectively. Performing renormalization consistent with the use of the NLO weak Hamiltonian 
and its renormalization-scheme independence |]15|, the result, as expected, is given by 
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3T 



(3.9) 



where AG has aheady been defined by (|2.4|) and the value of the parameter is determined from 
momentum distributions in the process to be M^. 

The contributions from other two penguinhke diagrams, Fig. |^(d), are lengthy because they 
involve the b d'~f*g* vertex calculation [|16| and will be given only in the final form, in expression 
( p. 171) , and in the Appendix. 



The next step to be performed is to express the hard scattering amplitudes, ( |3.4| , p.6| , p^ in the 
Fourier-transformed space of transverse momenta. The Fourier-transformed amplitudes read 



'-fact 



Cf^-^IMlr { [(1 + X3(l - r2))e ■ pg] /^^(Dg, A, 6i, 63) 



TT 



+ 



T 



nonfact 



+ 



T 



nonfact 



Xi(l - r^)e • pi - xie ■ p^) hBiDc, D^^^bs, 61)} 

(2xi -X2-X3- r^(x2 - X3))e ■ ps hoiDq^, Dq, ^2, &i)} 
«,(/i) CpN, 1 



hc{Dg^,DG,b2,bi) 



TT 



-^^32M|r { [(1 - 2xi - X2 + X3 + r\l - X2 - X3))e ■ p^ 



hciDq^, DG,b2,bi) + {1 - r^){xi - X2)e ■ pi /^^(-Dgg, D^, 62, &i)} 



(3.10) 



where 



hAiDa, A, 61, 63) = Ko ^-q'albil K, J -ql\b^ + 63I 5{b 



hB{DG,D,^MM) = Ko [^\J-qh\h\j Ko -q',i\bi + b^lj ^(62) , 
hc{Dg„ Dg, b2, 61) = Ko (xTqllhl) Ko (J^\h - hi) 6{h + h 



hoiD^,, Dg, 62, &i) = Ko ( V-^3l^2| ) Ko ( ^-ql\bi - 62I ) S{b, + 63) 



(3.11) 



and Ko is the modified Bessel function of order zero. 

For the scales /i appearing in formulas ( p.lO| ) and ( p.9|) we are going to take the largest mass 
scale in the particular diagram : 



tA = max (y\l-qG, \/-^, l/^i, , t 



B = max 



tc = max 



to = max 



(\/^.\/^'V&l,l/&2) , 



oop 



max J-qc, I/&1 



(3.12) 



which ensures the reliable perturbative calculations with the small cts-coupling. 

Performing trivial b integrals over S functions and performing angular integrations by using 
Graph's theorem: 



(3.13) 



/(x, 61, 62) = J d(j)Ko{x\bi ± 62I) 

= 27r[e(6i - b2)Ko{xb^)Io{xb2) + 6(62 - b^)Ko{xb2)Io{xb^)] 
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one can finally write the total amplitude of the B ^ n (p decay as 



V ^ q=u,c 

Ag = |3(C3 + C5 - ^(C7 + Cg)) + C4 + - ^ (cg + Cio)| {Tfact) 
+ (C4 — -^Ciq) {Tnonfact{C>4)) + (Ce — -^Cs) {Tnonf active)) 



2ap 



V 9 



■ C2{Tioop) 



Q ' 



(3.14) 



with the matrix elements 



{Tfact) = — ^UM^r J dxidxs J bidbib3db3^B{xi,bi)^l{x3,b3 



X 



TT 



-H{Dg,D,MM) (l + X3(l-r2))e.p3 



,-(SB(tA) + ^.(tA)) 



TT 



(T„o„/,,i(04)) = -^UMlr I [dx] I bidb, b2db2^B{xu Ki^3, b 



H{Dg, ^gi, 63, bi) - r^)e ■ pi - Xie ■ p3 

Kix2,b2) 



X 



as(t 



c) 



TT 
TT 



H{Dg^, Db, &2, bi) (1 - r^)(l - xi - X2)e ■ pi 



^(5ij(tc7) + 5^(tc) + 5^(tc))|63 



H{Dg.^, Dh, 62, &i) (2xi -X2~X3- r^(x2 - X3))e ■ P3 



'iSBito) + S^itn) + S^itD))\b,=b, 



{TnonfactiO^)) 



X 



3 

as{tc) 

TT 



f^Mlr [dx] bidbib2db2^B{xiM 



An 



K{x3M) 



H{Dq^, Db, 62, bi) (1 - 2xi -X2 + X3 + r^(l - X2 - X3))e ■ p3 



iSBitc) + S^itc) + S^{tc))\b, 



--bi 



+ ^^H{Dg„ Db, 62, &i) f(l - r^)(xi - X2)e ■ p. 



TT 



^-{SBito) + S^{tD) + S.j,(tD))\b-i=bi 



where the Fourier-transformed expressions for the propagators in Eqs. ( p.l5| ) have the general 
form 



(3.16) 



H{D,,D2,bi,b2) = Ko{^-D,bi)f{^-D2,bi,b2) , 

with the function / as defined by (|3.13|) . 

In ( |3.14| ) we have also included the one-loop contributions from the diagrams of Figs, ^(c) and 
^(d). The matrix elements receiving the contributions from the diagrams in Fig. |^(d) are 

{Tioop)q = — ^f^Mlr J dxidx3 J bidbi <l>B(a;i, 6i)$* (xa, 61) 

^-{SB{tloop) + ST,{tioop))\b:i=bi I 



X 



1 ] as{Uoop) 



Ml 



Kn 



TT 



(3.17) 
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The expressions for and are given explicitly in the Appendix. 

One should note that in the above expressions we have pulled out the normalization factor 
/^/2y6 (/^ is the (/)-meson decay constant) of the wave function and we denote the rest by 
in order to have the same prefactor in both factorizable and nonfactorizable contributions. 

From the expression for {Tfact) in Eq.( p.l5D it is easy to essentially recognize the factorization 
structure in which a matrix element of a four-quark operator factorizes in the product of two 
current matrix elements (0|(ss)y_yi|O) ■ {n' {p'i)\{dh)v -a\B~ {pi)) ~ /^e^ ■ {Ap^ + -Bpa), and the 0- 
meson wave function integrates out. The current matrix element {'K~{p^)\{dh)v~.A\B~{pi)) exactly 
describes the B ^ n transtion form factor at the momentum transfer p"^ = {pi — Ps)^ = In 
Sec. IV we use this form factor to select mesonic wave functions. 

The expressions for the Sudakov exponents S^, S^, and Sb in Eqs.( |3.l5| ) and ( |3.17] ) are given 
in Sec. IV. 



IV. MESONIC WAVE FUNCTIONS AND SUDAKOV FACTORS 



The calculation of the matrix elements requires the knowledge of the scalar meson wave func- 
tions $. The hadronic wave functions represent the most speculative part of the perturbative 
approach. They are of nonperturbative origin and should be a universal, process-independent 
quantity. 

However, even for the most theoretically and experimentally exploited hadron, namely, the 
pion there are contradictory conclusions about the specific form of $7r- Theoretical calculations 
performed by using the QCD sum-rule method and on the lattice cannot distinguish 



between the most promising forms of the pion wave function, the asymptotic one 



<l>^'(x) = 6x(l - a;) 



(4.1) 



and the Chernyak-Zhitnitsky (CZ) wave function [1^ 



$^^(x) =30a;(l-x)(l-2x)2 



2VQ 



(4.2) 



both normalized to satisfy the experimentally obtained value for the pion decay constant = 
0.133GeV. Besides, any comparison between theoretically calculable processes and existing exper- 
iments cannot provide an unambiguous determination between them 
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What we are going to regard as the CZ wave function throughout the paper is the CZ form in 

1211: 



which the evolution from the hadronic scale /xq ~ 0.5GeV to some scale /ii is included 



l + (5(l-2x)2-l) 



a 



50/81' 



"s(/io) / 



2VQ 



(4.3) 



and, in the modified perturbative approach, the scale /ii is taken to be 1/6 pT| . 

If we are going to retain the intristic b dependence of the wave functions $ in the expressions for 
the matrix elements, then we are faced with even more uncertainties coming from the ambiguity 
in the form of the wave function 6-part as well as in the values of some new parameters. 

The constituent quark model of the wave function associates some Gaussian exponential to the 
6-dependent part ||22|, so that 
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b) = 6x(l - x)47r exp - x)6V(4c 



2^6 



(4.4) 



and 



<l>f^(x,6,/ii) = 6x{l - x) 



50/81' 



An exp (^— x(l — x)6^/ (4a^^) 



2v^ 



(4.5) 



where the pion's transverse parameters aas and are fixed from the tt — > 77 process to be 
aas = OMQGeV-^ and acz = O.Q55GeV'\ respectively |23|. 

For a S-meson wave function there exist a few models p4[. We consider two forms that have 
been proved in the calculations of various nonleptonic B decays. The first one is 



C + 



ml 



+ 



X x{l — x) 



(4.6) 



whose Fourier transform gives 



^%\x,h) 



iV(i) 



hx'^il — x) 



J Mix + Cx{l - 



(4.7) 



with the approximation m^ ~ = 5.286*61^. Ki is the modified Bessel function of order one. 
Neglecting the h dependence leads to 



A^W x{l-xf 
167r2 Ml + C(l - x) 



(4i 



For the constants A^*^^^ and C we have used the fitted parameters A^*^^^ = 604.34GeT^^ and C 
—27.5GeV^, which have been proved in other calculations . 



Another model is the oscillatorlike wave function of Bauer, Stech, and Wirbel p7 



iV(2) 

2n 



x{l — x) exp 



2u 



-X exp 



(4.9) 



with the constants A^^^) = l5QMGeV and tu = OAGeV. Both wave functions, ( [4.7|) and 
normalized with = 200MeV. 



are 



The vector-meson wave functions are modeled in the QCD sum-rule calculations |jT^, p8 



Since the form of the 0-meson wave function is still questionable, we have decided to use the 
asymptotic form 



^4>{x) = 6x{l - x] 



2^6 



0.233^61^. 



(4.10) 



without including any b dependence. We believe that owing to the lack of better experimental 
data to which transverse parameters can be fixed, an unrealistic 6-dependent part may produce 
more questionable results than by neglecting it. 

In order to suppress the soft contributions in the hard scattering amplitudes ( |3.15| , |3.17|) , we 
have included the Sudakov factors They ensure that the hard scattering amplitude receives 
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contributions only from the exchange of hard gluons, suppressing the contributions of soft gluons 
from the large b region. The Sudakov suppression is comprised by the hadron wave function 
redefinition 



$^ $0(x2,fe2) exp{-S^{t)) , 

^ $,(X3, fes) exp(-S^(t)) . (4.11) 

The Sudakov exponentials exhibit the result of all-order resummation of double logs appearing 
from the overlap of coUinear and soft divergences |29|. In our case, 



Ssit) = s{xip^ ,bi,t) - l/po In 



[ln{l/{b,AQCD)), 



S,{t) = s{x,pt,b,,t) + sm-X2)pt),h,t) - l//3oln f . , 

\ln{l/{b2AQCD))J 

SAt) = s{x,p^M + sm - x,)p^)M - l//3oln (i^^l^l^^^^) , (4.12) 

where Pq = (33 — 2nj)/12 and nj = 4. For Aqcd we have used the value Aqcd = 0.2GeV 
throughout the paper. The last term in the above expressions accounts for the renormalization 
from the IR scale 1/6 to the some renormalization scale t, which we are going to take to be one of 
the scales from ( |3.12| ), depending on the diagram considered. 

The full expressions for the Sudakov functions s{xi,bi,t), together with the usual approxima- 
tions used in a numerical treatment, can be found in ||30|| . 

Note that we have also associated the Sudakov function s{xiPi ,bi,t) with the light antiquark 
of the B meson. The heavy b quark, having a finite mass, does not produce collinear divergences 
and its Sudakov function is zero. 

Use of the above mentioned diversity of the wave functions would certainly diminished the 
capability of perturbative calculations for giving reliable predictions for the B~ 7r~0 branching 
ratio and CP asymmetry, having in mind that the effects of the large reduction of the results owing 
to the intrinsic 6-dependence in the wave functions as well as the large difference in the predictions 
depending on the B and tt meson wave function employed, has already been observed in other 
perturbative calculations ||3l|. We have checked that this is also the case in the calculation of the 



B ^ TT(j) decay. 

Owing to the specific character of the B ^ ncj) decay governed by the b — > dss transition where 
the strange quark- antiquark pair has to form the final 0-meson state, we can assume that the 
B Ticf) process is determined predominantly by the B ^ tt transition at the energy = M|. 

Therefore, we can try to make a selection among the wave functions by comparing the results 
for the B ^ IT transition form factor obtained from the QCD sum rule and lattice calculations 
summarized in 

F_f ^^(0) = 0.25 - 0.35 , (4.13) 

with those estimated in our modified perturbative approach. 

The expression for the form factor in the perturbative approach has the form 



[t]) = -Y^B J dxidxs J bidbib^db^i^BixiMWAxzM 
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H{Dg, D,, bu h) [1 + xsv] e-(^^(^^) + ^-(^^)) 
^^H{DG,D,„bM -r/)]e-(^^(^^) + ^-(^^))| , (4.14) 



which can be easily recognized in the expression for the factorizable part of the B ^ ircf) decay, 
Eq. (p.l5| ). The parameter r] is the fraction of the energy of the vr meson and at the momentum 



transfer = or = M| we have 77 = lor?7 = l — M|/M^ = 1 — r^, respectively. 

Estimating the B —>■ n transition form factor at the momentum transfer = using different 
forms of the B and vr meson wave functions taken from above, we achieve predictions which are far 
from the values obtained in the QCD sum rule and lattice calculations ( [4. 13] ), except if we assume 
the oscillatorlike model for the B meson wave function $^ (x), ( [4.9| ), and the CZ type of the pion 
wave function ([4.3| ), both being intrinsic b-independent. Our predicted value for the i? — >■ tt form 
factor obtained with these wave functions is 

Fpe7/iO) = 0.282 . (4.15) 

Both, $^^(x) and $^^(x, /xi) are more end-point concentrated wave functions than their alter- 
native forms, $^''(x) ( [4.8| ) and $"*(x) ( [4.1| ), respectively. This indicates a need for the enhance- 



ment of the soft contributions in order to match the predictions (f4.13 ) for the i? — > vr form factor 
esimated by nonperturbative methods. 

Comparable calculations of the 5 — > vr form factor in the modified perturbative approach 
have also been performed in |^ and, similarly, the results obtained have exibited strong 



dependence on the mesonic wave functions used, confirming that the wave functions represent the 
weakest point in the calculation of 5-meson decays in the perturbative approach. 

V. NUMERICAL RESULTS AND DISCUSSIONS 

Now we are going to discuss the branching ratio BR and the CP asymmetry in the B^ n^cf) 
decay numerically. 

The decay rate is given as 

where \^'^{Mb, M^, 0) = Ml{l-r^) and the total amplitude M is given by ( p.l4| ) . CP asymmetry 
in terms of the Au and Ac amplitudes, Eq.( |3.14| ), reads 

^ -2VJm{y^)Im{AuAl) 

'''''' (IKH + IKPIAP + 2VcRe{y^)Re{AuAl)) ' ^ ' ' 

The products of the CKM matrix elements may be written in the Wolfenstein parametrization 

as 

K = = A\\l - AV2)(p - iTi) ^ A\\-p -tfj), 

Vc = V:,Vcb = -AX' . (5.3) 

We use the following values of the parameters p and rj: 



15 



p = 0.16, r/ = 0.33, (5.4) 
which correspond to their central values obtained by the unitarity fit Since recent measure- 



ments disfavor the negative values for the p parameter , the CP asymmetries will be presented 
in figures by taking p in the range 



< p< 0.25 



(5.5) 



and using the central value for the rj parameter from (|5.4| ) , the minimum and the maximum allowed 
values 



V 



0.27 , 
0.38 , 



(5.6) 



respectively. The other Wolfenstein CKM parameters used are A = 0.823 and A = 0.2196. 



Following Ref . , we are going to take the constituent quark masses in the loop expressions, 
Eqs.( |2.4| ) and ( 3.17 ), with particular values = 0.2GeV and rric = 1.5GeV. 

For the Wilson scale-independent coefficients at the renormalization scale /i = ruf, = 4.8GeV 
(asiMz) = 0.118, aem(M^) = 1/128) we take [0 



ci = -0.324 
C3 = 0.017 , 
C5 = 0.011 , 



C7 



-1.05-10" 



C2 = 1.15 , 
C4 = -0.038 
C6 = -0.047 

cs = - 



-3.84 ■ 10"^ 



Cg = -0.0101 



cio = 1.96 • 10" 



(5.7) 



One can note from expression ( |5.2|) that some absorptive part in the amplitude is essential for 
nonvanishing CP asymmetry. 

As in the BSW factorization approach, the necessary absorptive part comes from the cut in 
the penguinlike diagrams in Fig.^(c), residing in the term AG[m'^, M|, /x^). From expression ( p. 41 ) 
it is easy to see that the absorptive part is developed for the virtual photon momentum q such 
that > 4mg, q = u,c. Owing to the specific momentum distributions in the B~ 7r~0 process, 
the imaginary part emerges only from the diagram with a u quark in the loop. In numerical 
calculations we use the approximation of (|2.4| ), 



AGappirrig, 



- + - + (l + -)R(z)-\n^ 
3 z z p^ 



(5. 



where by defining a = J\l — A/z], we have 



R{z) 



—an + 2a arctan(a) 
1 — a 

tan + a In ■ 



1 + a 



z = {M^/m^f 



(5.9) 



In addition, in the perturbative approach there are absorptive parts connected with the cuts in 
the propagators of virtual partons in each of the diagrams in Fig. |^. The expression for if, ( p. 16 ), 
can develop the imaginary part for some of the values of the fractions xi in the integration for 
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TABLE I. Branching ratios for the Ti^ip decay calculated for different penguin contributions, 

by taking only the factorizable parts or the complete expression into account. The first column contains 
the predictions obtained in the BSW factorization approach by using (g^) = M| (see text). Columns I 
and II give predictions calculated in the modified perturbative approach by employing the meson wave 



functions ^%\x), Eq.(|l|) and Eq.Q, and by using the Wilson coefficients Ck{mh) (K 



in column I and cf\t) ( |5.16| ) in column II. The CKM parameters used are p = 0.16 and r] = 0.33. 



Penguin contributions 




BR 




BSW 


I 


II 


QCD-factorizable 


0.20 • 10"^° 


0.14-10"^° 


1.06 • 10^1° 


QCD-aU 




2.51 • 10-1° 


0.73 • 10-10 


QCD+QED-factorizable 


0.34 • 10-*^ 


0.38 • 10-s 


0.89 • 10-8 


QCD+QED-all 




0.44 • 10-s 


0.85 • 10-8 



which the denominators of the gluon or quark propagators under the square root become negative 



see Eqs.(p.5|, |3l8|)). In this case we take 



Koityb) = '^H^^\yh) (5.10) 

and 

f{iyMM) = y [e(6i - h)H^^\yh,)Uyh2) + 6(62 - h)Hi^\yh2)Myhi)\ . (5.11) 

('2') 

Having selected the B meson wave function $y (x), ( |4.9|) , in the preceding section: 

^f{x) = ^^%\x. 0) = ^^x{l-x) exp (-^a;^) , (5.12) 

and the $^^(x, 1/6) ( ^.3| ) for the vr meson wave function, we can now continue along the lines 
developed in the preceding sections and give reliable predictions for the n^cf) branching ratio 

and CP asymmetry in the modified perturbative approach, using the NLO weak Hamiltonian. 

The results are presented in Tables I and II, together with the predictions estimated in the BSW 
factorization approach, both being calculated with the preferred values of the CKM parameters, 
p = 0.16 and rj = 0.33. 

Calculations of the n^cf) branching ratio and asymmetry in the BSW factorization 

approach have been performed by many authors [Q- 0. In order to be able to clearly assign 
the role of nonfactorizable contributions in the decay, we have recalculated the BSW factorization 
predictions using our values of the Wilson coefficients Cj ( |5.7| ) and the CKM parameters p and rj. 

The decay amplitude in the BSW approach can be directly compared with the complete ex- 
pression for the B -Kcf) amplitude given by Eq.( ^.14 ) by neglecting nonfactorizable parts and 



numerically suppressed contributions emerging from the diagrams in Fig. 3(d). The matrix ele- 
ment in the strict factorization approach is propotional to the F^^'^ (M^) form factor which we 
calculate in the single-pole approximation as 
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TABLE II. CP asymmetries for the B~ — > 7r~(/> decay calculated for the QCD and QED penguin 
contributions together, by taking only the factorizable parts or the complete expression into account. The 
first column contains the predictions obtained in the BSW factorization approach by using (g^) = M| (see 
text) . Columns I and II give predictions calculated in the modified perturbative approach by employing the 



meson wave functions ^>2^(x), Eq.(|5l2D and 1/6), Eq.(|^, and by using the Wilson coefficients 



Ck{fnb) (^3) in column I and cf\t) ( |5.16| ) in column II. The CKM parameters used are p = 0.16 and 
rj = 0.33. 



Penguin contributions 




acp/lO-^ 






BSW 


I 


II 


QCD+QED-factorizable 


-1.9 


14.6 


15.4 


QCD+QED-all 




16.1 


16.3 



^BSW^M^) - ^ _ M2/M|,(l-) ' ^^-^^^ 

where M^^(l^) = b.?>2GeV and F^^{0) = 0.33 It is worth mentioning that the prediction 
for the B ^ IT form factor estimated in the perturbative approach, Eq. ([4.15|) , is somewhat smaller 
than the F|^(0) value. 



All results estimated in the factorization approach are obtained by taking the virtual photon 
momentum squared equal to = M|. We have already stated that, in general, information about 
the g^-value is lost by factorizing hadronic matrix elements, except in the strict factorization, 
when nonfactorizable and/or strong final-state interactions are neglected. Therefore, considering 
possible nonfactorizable or final-state corrections, after the factorization procedure, the is usu- 
ally considered as a free parameter, whose average value is constrained by some simple, general 
kinematical reasons to be ||3^ 



ml/A < {q^) < ml/2, (5.14) 

and usually assumed to be valid for all nonleptonic heavy-to-light transitions. 

The dependence of the branching ratio and CP asymmetry on the (g^) as a function of the p 
CKM-parameter is shown in Figs. ^. The branching ratio appears to be practically independent 
of the value of the (g^). Such behavior is due to the cancellation which occurs between the 
Wilson coefficients multiplying the one-loop QED penguinlike matrix element, ( 3.14 ), because of 
the relation 

3ci(mb) ~ -C2{mb) . (5.15) 

On the contrary, CP asymmetry exibits a large reduction of up to 70% if higher (g^)-values are 
taken. 

The results in column II in Tables I and II are obtained by calculation inspired by the papers 
of Li and collaborators , , in which the fi scale-setting ambiguity of Wilson coefficients is 



moderated by applying the three-scale factorization theorem. Their theorem keeps trace of all three 
scales characterizing the nonleptonic weak decay, the W-boson mass Mw, the typical scale t of the 
process, and the hadronic scale ~ Aqqd, and proves for the leading-order weak Hamiltonian that 
Wilson coefficients should be taken at the scale t, a typical scale in a particular decay. The matrix 
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0.0 0.1 0.2 0.0 0.1 0.2 
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FIG. 4. Branching ratio and CP asymmetry in the — > 7r~4> decay calculated in the BSW 

factorization approach as a function of the CKM parameter p and for the central value of ?7 = 0.33. The 

solid, long dashed, and dot-dashed lines denote predictions obtained by taking (g^) = M|, {q'^) = 

and {q'^) = respectively (see text). 



elements of the operators Ok and Wilson coefficients are then both calculated at the same scale. 
The scale is determined by the dynamics of the process, contrary to the arbitrary renormalization 
scale fj, taken to be a constant, nib, for the Ck coefficients, ( |5.7| ). 

Under the assumption that the three-scale factorization theorem is also valid for the NLO weak 
Hamiltonian ( p.lD we have taken the explicit form of the Wilson coefficients, calculated directly 
at Mw and then rescaled to some lower scale t [H]: 



JO) 



C5 



7.(0) 



.(0) 
-10 



1 + 0{as{t)) + 0{ae 



as{t) 



247r 


as{t) 




Stt 






247r 


asit) 




8tt 






4 


Gtt 




0, 




Q^em 


4 


Gtt 




0, 





Eo{xt) + 2 lo? 



10 



+ 



1 



2 

i^o(x0 + 3log^ 



w 



10 



^o(a;t) + 77 log 



3 M,2 



10 

e 10 

9" 



4, 

T^log 



AC^ixt) + Do{xt) + - log 



9 



w 



Gtt sin 6vk 



[25o(xi)+Co(xi)] 



20 
27 

20 
27 



+ 



1 



Stt sin 014/ 



[hB^{xt) - 2Co{xt)] 



(5.16) 



where Xt = mf/M^r. The functions Bq, Cq, -Dq, and £"0 are the Inami-Lim functions . 
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These coefficients are only an approximation of the Wilson coefficients Cfc(/i = rrib), 
obtained by performing the renormalization-group analysis [|T1| and used throughout the paper, 
but we hope that possible uncertainties involved in the calculation by using the coefficients in (^.16|) 
are covered within the accuracy of our model. In the numerical estimates we have also taken into 
account 0{as) corrections in cf^ and 0^2^ in order to have the proper O(asaem) calculation. 

By taking the Wilson coefficients cj^^ at one of the scales ( p.l2| ), depending on the diagram 
involved as a contribution of the operator Ok, we obtain the results given in column II in Tables I 
and II. The results are estimated again with the selected wave functions ^^^\x) and ^^^{x, 1/6). 

Let us now discuss the results from Tables and emphasize their general characteristics. 

One can note that the -K^cf) process is clearly dominated by the EW penguin contribu- 

tions, in both the factorization and the perturbative approaches and the predicted branching ratio 
for the B~ — >■ 7r~(/) decay is of order 0(10"'^). 

It is obvious that EW nonfactorizable contributions are small, being directly proportional to the 
small Wilson coefficients cg and ciq. For nonfactorizable contributions of QCD penguin operators 
there is no such apparent reason, because the Wilson coefficients multiplying the operators C4 and 
Ce, are in absolute magnitude even larger than the coefficient Cg, which dominates the 7T~(j) 
decay ( |5.7|) . The influence of the QCD nonfactorizable contributions is noticeable, specially in 
the perturbative results based on the Wilson coefficients taken from ( ^.7|) and represented in 
column I. For this case, by comparing the third and fourth rows in Table I, we can see that 
nonfactorizabile corrections can account for some 14% of the final result. However, after taking 
the Wilson coefficients in the convolution with the hadronic matrix elements at the same scale, as 
it is done by obtaining the results in column II, nonfactorizable contributions become negative, 
and small, and can be considered negligible, lowering the final result by some 4%. Negligible 
nonfactorizable corrections in this model indicate that, by a suitably chosen scale which truly 
makes the product of the Wilson coefficients and the hadronic matrix elements scale independent, 
it is possible to account for the almost strict factorization in the B^ ix^cj) decay, which would 
be naively expected by the "color transparency argument" p|. 

Further general behavior of the results from column II can be summarized in the statement that 
the branching ratios calculated using the Wilson coefficients c^^"* (t) are enlarged by some factor two 
in comparison with the estimations obtained by using Cfc(yU = mi,). In addition, CP asymmetries 
are predicted to be about 16%, similarly as in the model in column I, where Cfc(m{,) are used. 

The predicted asymmetries are much larger than those obtained from the BSW factorization 
model, and they are of an opposite sign. The reason for such an enlargement of CP asymmetry are 
absorptive contributions due to the on-shell effects in the propagators of virtual partons appearing 
in the perturbative calculation. They are also present in the factorizable amplitude {Tfact) which 
multiplies the penguin loop. Neglecting of the imaginary parts coming from the on-shell effects in 
the propagators would give predictions for CP asymmetry comparable with that obtained in the 
factorization approach. 

Figures ^ show the impact of different choices of CKM parameters on our predictions for 
the branching ratios and CP asymmetry in the B^ 7i^(f) decay. The predicted asymmetries 
calculated by using the Wilson coefficients Ck{fnb) and c^k\'t) are almost the same. Therefore, we 
show explicitly only the CP asymmetry obtained by using c^k\t) coefficients. One can note that 
the predicited CP asymmetry can be enlarged up to 22%. 
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FIG. 5. Branching ratio and CP asymmetry in the TT^(j) decay calculated in the modified 

perturbative approach as a function of the CKM parameter p. The solid, long dashed, and dot-dashed 
lines correspond to the values of the CKM parameter rj = 0.33, fj = 0.27, and r] = 0.38, respectively. 
Predictions obtained by using the Wilson coefficients ci:{rrn,) ( |5.7D and c^j^\t) ( |5.16| ) are denoted by labels 
I and II, respectively. 



VI. CONCLUSIONS 



In this paper we have calculated the branching ratio and CP asymmetry of the penguin- induced 
7T^(f) decay in the modified perturbative approach by applying the NLO effective weak 
Hamiltonian. Working in the framework of the modified perturbative approach we have included 
the transverse momentum dependence and the Sudakov form factors. The modified perturbative 
approach also enables us to calculate nonfactorizable contributions. 

We have used the B ^ n transition form factor to select mesonic wave functions by com- 
paring our result with the predictions estimated in the QCD sum rule and lattice calculations. 
The comparable prediction has been obtained only for the intrinsic b-independent, more end- 
point concentrated wave functions for both B and vr mesons, $g)(x) and ^^^(x, 1/6) f^), 
respectively. 

Using the NLO weak Hamiltonian and the selected wave functions we have first worked with the 
renormalization scheme-independent coefficients and have been able to calculate the EW penguin 
contributions properly, proving their dominance in the B~ 7r^0 decay. 

In addition, we have examined the assumption of taking the Wilson coefficients to be convo- 
lution functions in the starting factorization formula ( |2.5|) instead of taking them as constants at 
some arbitrary scale /x. The Wilson coefficients then enter into the factorization formula in the 
convolution with the matrix elements at the same scale t, typical of the process and that resolves 
the problem of different renormalization scales for the short-distance part (Wilson coefficients) and 
the long-distance part (matrix elements of the four-quark operators) in the amplitude of the weak 
Hamiltonian. Estimations based on this assumption have produced the branching ratios about 
factor two larger than those calculated with the conventional Wilson coefficients. 

Besides, if the Wilson coefficients are considered to be functions of the scale, the same one which 
appears in the hadronic matrix elements, then the nonfactorizable QCD penguin contributions 
appear to be negligible, as is the case with the obviously very small nonfactorizable contributions 
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of the EW penguin operators. 

Therefore, our results for the branching ratio appear to be in agreement with previous calcula- 
tions performed in the BSW factorization approach, predicting the branching ratio to be of order 
O{10~^), dominated by the EW penguins. On the other side, the predicted CP asymmetry differs 
a lot from that estimated in the BSW factorization approach, being as large as 16% and having 
an opposite sign for the preferred values of the CKM parameters p = 0.16 and rj = 0.33. The large 
CP asymmetry estimated in the perturbative approach is the result of large on-shell effects of the 
virtual propagators involved in the calculation. 

The strong reduction of the results obtained with the intrinsic b-dependence of the wave func- 
tions indicates that mesonic wave functions still needs further investigations. Presently, i?-meson 
wave functions suffer from uncertaintes involved in the models from which they are derived as 
well as from uncertaintes coming from the fit to experimental data, and ask for a more refined 
treatment in their derivation. 

Provided that the 5-meson wave function could be better determined, the formalism of this 
paper may be successfully applied to similar penguin-induced decays, of which K~(f) and 

B~ —>■ uj(f) are particularly interesting owing to recent experimental measurements and their role 
in the determination of the values of some CKM matrix elements. These topics will be the subject 
of our future investigations. 
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APPENDIX A: 

The results of the calculation of diagrams E and F, shown in Fig. ^(d), are given here explicitly. 
In the calculation we have neglected the transverse momenta in the loops. The general expression 
for the contribution of the diagrams can be written as 

T^=C' + II + H; + GI , t = E,F. (Al) 

The results are presented as integrals over Feynman parameters and the particular contributions 
are found to be 

= (-)^[2(1 - r2)(l - 2xi)e ■ pi + (1 + - 4xi + 2x3(1 - r^))e ■ ps] , 
6 

du / dv(—) — — (u + v(l — Xi))(l — u — v(l — Xi))(l — u — V + (1 — v)(—2xi + X3) 

Jo Mr. 

hi 



+ r2(l -u-v~ 0:3(1 - v))) [(1 - r^)e ■ pi ~ e ■ ps] , (A2) 

rl rl—u Tfl^ 

= du dv^\{l-r'^){u + v{l-xi))e-pi 
Jo Jo Me 

,2 



+ (1 -2u-2v~ 2xi(l -v) + x^il - r'))e ■ p^] , (A3) 



rl rl-u m'^ 

Gf = du dv\n^\{l~r^){l-2u-2v{l-xi))e-pi 

^ Jo Jo jJL^ 



-(1 -u-v - 2xi{l -v) + X3(l -v) + r^(l - u-v - x^{l - v)))e -ps] , (A4) 
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and 



;[(l-r2)(l -2a;i)e-pi + 2(1 + r^- 4x1 + 2x3(1 -r2))e- Pa] , 



du 







, Ml 
dv — — 
Mi 



X |(1 — r^){—uxi + ^(1 — Xi))[(f (1 — Xi) + Xi{l — u)){{l — f )(1 — xi) + ■ua;i)e ■ pi 
+ ((w(l - t;) + a;i(l - M - 3f + 2f (m + v))){l - - 2((1 -u){l-u- 2v) + v^)xia;3 
+ {2{u -v) + {u + vf -{{l-u){l-u- 2v) + vy^)xl)e ■ P3]} , (A5) 



Gi 



du 



l-u 



dv 







Ml 



-[(1 - r^)xit ■ pi - {{v{l - X3) + a;3(l - u)){l - r^) + r^)e ■ ps] , (A6) 



— du 
Jo Jo 

+ [2{xi -X3){l-u-v)- r^{l -2v- 2x3(1 -u- v))]e ■ ps} . 



dv In — ^ 
//2 



^{1 — r^){v — Xi{u + v))e ■ pi 



(A7) 



The functions M^ and Mp depend on the quark mass in the loop and are given by 

Ml = ml- M|[m(1 - uy + v{l -v){l- xi){-{xi - x^) + r^{l - X3)) 



-uv{-{xi - X3) + r^(2 - xi- X3))] 



(A8) 



and 



Ml ^rrl-Ml 



u{l - u)xi{xi - X3 + r^X3) + f (1 -v){l- xi){-{xi - X3) + r^(l - X3)) 



- uv{{xi - X'i){l - 2xi) + r^(xi + X3 - 2x1X3)) 



(A9) 
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